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Abstract The flow of superfluid helium at very low temperatures around an oscillating 
microsphere is known to be unstable slightly above the critical velocity. The flow pattern 
switches intermittently between potential flow and turbulence. From time series recorded 
at constant temperature and driving force the statistical properties of the switching phe- 
nomenon are discussed. Based on our recent understanding of the critical velocity v c of 
5-H 



oscillatory superflows, i.e., v c ~ \/KCO, where K is the circulation quantum and CO is the 
oscillation frequency, the analysis is being refined now. From the exponential distribution 
of the lifetimes of the turbulent phases quantitative information on the vortex density L can 
be inferred such as the distribution and the width of the fluctuations of L. The phases of 
potential flow show a Rayleigh distribution of the excess oscillation amplitude above the 
amplitude at turbulent flow. The rms value is found to scale as a>~ 3 ' 2 . 
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1 Introduction 



> 

Os. 

^^ , Laminar and turbulent flow of superfluid He around an oscillating microsphere (radius r = 

0.12 mm) has been studied in detail in recent years Q]|2). At small oscillation amplitudes 
the drag force is linear in velocity amplitude (Stokes regime) and is determined by the nor- 

r — . mal fluid component. At temperatures below ca. 0.5 K the drag is given solely by ballistic 

phonon scattering and vanishes as T 4 . Above a critical velocity amplitude v c a transition to 
a large and nonlinear drag force is observed that scales as (v 2 — v 2 ), see Fig.l. (The small 
remaining phonon drag can subtracted from the total drag force and will not be considered 
here any more.) This transition is interpreted as the onset of turbulent flow of the superfluid. 
For velocities just above v c an instability of the flow is observed at temperatures below ca. 
0.5 K, namely an intermittent switching of the flow pattern between turbulent and laminar 
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Fig. 1 Velocity amplitude v of the oscillating sphere (frequency 1 14 Hz) as a function of the driving force 
F at 300 mK. At small drives the flow is laminar and the linear behavior is determined by ballistic phonon 
scattering. The shaded area above the critical velocity v c ~ 20 mm/s indicates the unstable regime where the 
flow switches intermittently between potential flow and turbulence. At larger driving forces we observe stable 
nonlinear turbulent drag where F(v) scales as (v 2 — v 2 ). 



phases (2). 

In Fig. 2 time series recorded at constant temperatures and at three different driving 
forces are shown. During a turbulent phase the drag is given by the turbulent drag force and, 
hence, the velocity amplitude is low. It is obvious that the lifetime of the turbulent phases 
grows very rapidly with the driving force. When turbulence breaks down the velocity ampli- 
tude begins to grow due to the much smaller phonon drag. Because the velocity amplitudes 
are above v c these phases of laminar flow (better: potential flow) usually break down and a 
new turbulent phase is observed. In some cases, in particular at low drive, metastable po- 
tential flow with very long lifetimes is observed occasionally when the maximum amplitude 
corresponding to phonon drag is reached (see upper trace in Fig. 2), as will be discussed later. 



The statistical properties of this switching phenomenon have been analyzed in some 
detail already in earlier work J2][3]. Since then, our understanding of the transition to turbu- 
lence of oscillatory flow at a critical velocity v c has been greatly improved by identifying 
a universal critical velocity that does not depend on any geometry, namely v c ~ \/k- co, 
where K is the circulation quantum and co is the oscillation frequency |4.5]. This result can 
be derived from a general argument based on the "superfluid Reynolds number" and in more 
detail from the dynamical behaviour of the turbulent vortex tangle as recently discussed by 



Statistics of Intermittent Switching 



30 
28 
(/> 26 
E 24 
E 22 
^20 











time (100 s / division) ► 

Fig. 2 Three time series of the velocity amplitude at 300 mK and at three different driving forces (in pN): 47, 
55, 75 (from top to bottom). The low level v, corresponds to turbulent flow while the increase of the velocity 
amplitude occurs during a laminar phase, occasionally reaching the maximum value v max given by the linear 
drag. With increasing drive the lifetimes of the turbulent phases grow rapidly. 



Kopnin I4ll6l . It is the purpose of this work to present now a more refined analysis of the 
statistics of this switching phenomenon. In the following Section 2 the tools necessary for 
the statistical analysis will be outlined, in Section 3 the analysis of the statistics of the turbu- 
lent phases will be discussed. It will be shown that quantitative information on the turbulent 
vortex density L and its fluctuations can be inferred from the data. Moreover, it will be em- 
phasized again that superfluid turbulence is inherently unstable and only appears to be stable 
due to the finite time of experimental observation. In Section 4 the statistical properties of 
the phases of potential flow will be presented. Finally, the observation of metastable poten- 
tial flow above v r will be discussed. Section 5 concludes this work. 



2 Statistical Tools 



As presented earlier |2] the reliability function R(x) is the relevant quantity of our analysis. It 
is related to the cumulative distribution function F(x) which is the probability that a random 
variable has a value not larger than x: 



R( x ) = 1 - F(x) 



The failure rate A (x) is defined by 



dR 
dx 



-A(x)-R(x) 
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Fig. 3 (From Ref.3) Statistical analysis of the lifetimes of the turbulent phases, a) The exponential decay of 
the reliability function R(t ) gives the mean lifetime, b) The mean lifetimes measured at different temperatures 
and driving forces as a function of the turbulent drag force. The solid line is a fit of an exponential increase 
with the square of the driving force, see Eq.3. 



which gives 



A(x)- 



dlnR 
dx ' 



(1) 



We note that A is the conditional probability for an event to occur in the interval (x,x + 
dx) provided it has not occurred until x. A is a very important quantity because it contains the 
physics. For a constant failure rate, e.g., we obtain the exponential decay of R known from 
radioactive decay. For a so called Rayleigh distribution R(x) = exp(— x 2 ) we have A = 2x. 
With these tools we now evaluate first the statistics of the turbulent phases. 



3 Statistics of the Turbulent Phases 



In Fig. 3a the distribution R(t ) of the lifetimes of the turbulent phases is shown for a particular 
drive and temperature. From the straight line behavior we obtain an exponential decay 



R(t) = exp(-Vf) 



(2) 



with 1/v =< t > being the mean lifetime. Thus the failure rate A = v is constant, i.e., 
there is no history dependence in the probability for the decay of a turbulent phase. The 
values for the mean lifetime are independent of temperature and are found to collapse to a 
universal drive dependence if the strongly temperature dependent phonon drag Xv, at the 
velocity v, is subtracted from the external driving force. (A is the linear drag coefficient 
due to phonon scattering.) This is reasonable because the phonon scattering is probably not 
affected by the turbulence. Obviously, the mean lifetime <t > depends only on the turbulent 
drag force F nir t, = F — Xv t , see Fig. 3b. The rapid increase of the lifetime can be described 
by an exponential growth 
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Fig. 4 (From Ref.3) Sketch of a stationary normal process. The statistical fluctuations around a stationary 
mean are assumed to have a normal distribution. The down crossings of a given level below the mean are 
indicated. 



with t = 0.5 s and F\ = 18.3 pN. 



<t>=t exp(F lurb /Fi) 



(3) 



From Fig.l we find the turbulent drag force being given by F tur \, = y(v 2 — v 2 ), where 
7= coPs{T)7tr 2 /2 and C£> = 0.36 for our sphere 1T||2). F° r simplicity, from here on we shall 
omit the subscript "turb", i.e., assuming that the linear drag force is always subtracted from 
the external driving force (but keeping in mind that, strictly speaking, a correction factor 
of 8/371 as 0.85 is needed when comparing the driving force with the nonlinear turbulent 
drag force J7)): F(v) = y(v 2 — v 2 ). Writing yv 2 = Fq gives F(v)/Fq = (v 2 — v 2 )/v 2 and from 
L = (v/fc) 2 , where the amplitude of the vortex density L is the steady state solution of the 
Vinen equation_J|6], we have: 

£-^=£. m 

Fq L c L c 



where L c = (v c /k) 2 ~ (^/kco/k) 2 ~ co/k. We note that L c is the lowest vortex density for 
turbulence to exist and is independent of the drive. Any fluctuation of L that drops below 
L c will cause a breakdown of turbulence. For a typical drive F = 30 pN (see Fig. 3b) and a 
typical F = 450 pN (at 120 Hz) we have 



AL 



= — = 6.7-10 



-2 



AL 



= 6.3 10 



-2 



(5) 



F L F + F 

This means, that in this case a fluctuation larger than 6.3% of L will cross the level of L, 
downwards and, therefore, cause a breakdown of turbulence. How often does that happen? 



In the following a model is being presented that describes both the exponential distri- 
bution of R(t), see Eq.2, and the dependence of the mean lifetime <t>, see Eq.3. Let us 
consider a stationary normal process as sketched in Fig.4. The properties of level crossing 
by a random process are mathematically very well investigated [8|. First of all, the down 
crossings constitute a Poisson process on the time axis and, therefore, the time intervals are 



1 Here we apply the Vinen equation to oscillatory flow provided that the relaxation time of L is short 
compared to the period of the oscillations, for details, see (4). 
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exponentially distributed. We postulate that the lifetime of a turbulent phase is given by the 
first down crossing. Consequently, the times of the first down crossings are also exponen- 
tially distributed when the system is restarted many times with a new turbulent phase. This 
explains the exponential distribution of R(t ), see Eq.2. It is also plausible, that an increase of 
the drive will shift the level L further up and, hence, a larger negative fluctuation is needed 
for a first down crossing below L c and, consequently, the lifetimes of the turbulent phases 
will grow. For a quantitative analysis we apply Rice's formula [8,9| for the average number 

< N > of down crossings per unit time of a level C below the mean of a stationary normal 
process <§(f) with zero mean and variance a 2 , namely, 

<A.>^ v ^(0)exp-(-|-) 2 , (6) 

where r"(Q) is the second derivative of the autocorrelation function of the random process 
at the origin that is related to the second spectral moment by 

-r"(0) = / w 2 S{a)dco, (7) 

Jo 

with S(co) being the spectral density function of the process. The exponential dependence of 

< N > on the level C reflects the assumed normal probability density function of the under- 
lying process /(£) = (l/\/2Jr<7)exp(— E, 2 /2a 2 ). From the measurements of the average 
time < t > of the first downcrossing we have the condition 



and therefore 



<N> ■ <t>= 1 



In ( C \ 2 (F^ : 



Writing 



<t> =^ = 7^m ew \v^) = ' oexp UJ ■ (8) 

F F F AL F Q AL F x 

— = — ■ — = — ■ — = — , where L\ —L c - — , 

F\ F F x L c Fi L\ F 

we can compare the exponents: C = AL and \/2 a = L\ and finally we obtain for the width 
a of the Gaussian fluctuations 

_ U_ _ L C F X <J L\ _ Fi 

a ~ V2~ V2F an L ~ ^2L~ V2{F + F )' 

For a drive F = 30 pN and for F Q = 450 pN and F x = 18 pN we find a/L = 2.7 ■ 10~ 2 . In 
Fig. 5 this numerical example is sketched graphically. 

The prefactor to = 0.5 s in Eq.8 is interpreted as twice the average time between two 
consecutive crossings of level L |9|, i.e., there are 4 crossings per second. Hence, the aver- 
age frequency of the fluctuations is 1.5 Hz, which is surprisingly low. 

We now have achieved a quantitative interpretation of the statistical properties of the 
turbulent phases that no longer requires qualitative dimensional considerations as earlier 
|3]. Our model of a stationary normal process on which the Rice formula is based, allows a 
detailed analysis of our experimental results. An interesting implication is that although the 
lifetimes of the turbulent phases grow rapidly with increasing drive they never diverge. That 
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turbulence breaks down 



Fig. 5 (Color online) At a drive of F = 30 pN the vortex density L is 6.3% above the critical value L c and 
the relative width of L is 2.7%. All fluctuations crossing the level L c will cause a breakdown of the turbu- 
lent state. With increasing drive the level L increases linearly while L c and the width o~ remain unchanged. 
Consequently, the lifetimes of the turbulent state grow rapidly. 



means that turbulence is inherently unstable. When the finite time of experimental observa- 
tion is shorter than the lifetime a turbulent phase, turbulence only appears to be stable. An 
analogous situation has recently been observed in classical pipe flow where turbulence is 
shown to have lifetimes that grow very rapidly with the Reynolds number without diverging 
at some critical value 1 101. 



4 Statistics of the Phases of Potential Flow 



In Fig. 6 the statical properties of the laminar phases are displayed. The reliability function 
R is shown in Fig. 6a as a function of the excess velocity amplitude Av = v — v, above the 
velocity v, at turbulent drag at a particular temperature and drive. The data can be fitted by 
a Rayleigh distribution R(Av) — exp(— (Av/v w ) 2 ) and the fitting parameter v w is shown in 
Fig. 6b to be independent of temperature and driving force (phonon drag being subtracted 
again). We find, however, that v M . has some frequency dependence. For further discussion 
we are changing variables now by considering oscillation amplitudes Aa = Av/co rather 
than Av, i.e., we set a w = v w /a>. It is a property of the Rayleigh distribution that the fitting 
parameters v w and a w are the rms values of the reliability functions R. Although the range 
of frequencies of our data is rather limited, a power law fit to the data in Fig. 7 seems to 
indicate an 0)~ 3 ' 2 dependence of a w . The rms values a w drop from 9 /im at 114 Hz to 3.4 
/im at 212 Hz. This means that at high frequencies the laminar phases are breaking down at 
lower excess amplitudes Aa than at low frequencies. 

Discussing these results the following questions have to be considered: 1 . Why do we 
observe a Rayleigh distribution of the excess amplitudes Aa] 2. Why has the rms value a w 
this frequency dependence? 

The answer to the first question is made plausible by considering the failure rate A of 
the Rayleigh distribution A = —d XnRI dAa = 2Aa/ ' a\, that is proportional to Aa. As the 
sphere increases its oscillation amplitude during a laminar phase it will eventually encounter 
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Fig. 6 Statistical properties of the laminar phases, a) The reliabilty function R(Av) = exp(— (Av/v„) 2 ) is a 
Rayleigh distribution, b) The fitting parameter v w is independent of temperature and drive. 
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Fig. 7 (Color online) The rms amplitudes a w for six different oscillation frequencies. Each data point is 
obtained from the average of several individual time series recorded at different temperatures and driving 
forces. The error bars are calculated from the standard deviation. A power law fit that takes into account the 
error bars indicates a slope of -1.43 ± 0.14. 



remanent vortices that were trapped when the preceeding turbulent phase has broken down. 
The collision of the sphere with these vortices will give rise to a breakdown of the potential 
flow. (A similar observation has been made recently by the Osaka group with a vibrating 
wire colliding with vortex rings causing the flow to become turbulent lUll ). The probability 
for this event to occur is proportional to Aa. This situation is analogous to the hysteresis 
observed at higher temperatures where during the initial up-sweep of the driving force (after 
the cell has been filled with helium) the critical velocity v c can be exceeded considerably 
before the transition to turbulence finally occurs because the large normal fluid drag at tem- 
peratures above ca. 0.5 K limits the excess amplitude and hence the failure rate. A substantial 
increase of the drive is needed for the laminar phase to break down. Therefore, the switch- 
ing phenomenon is replaced by the hysteresis (2][5). It also fits nicely into this picture that 
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metastable potential flow can be observed when the maximum oscillation amplitude due to 
phonon drag is reached. In this case the spacing of the remanent vortices is larger than the 
maximum amplitude of the sphere, no vortices can be reached and, hence, potential flow is 
metastable. Only because of vorticity generated by the absorption of particles from natural 
radioactive decay the mean lifetime of these metastable states of superfluid potential flow is 
limited to 25 minutes (2]|7). 

At present, the second question can be discussed only on a qualitative level. The fun- 
damental length scale of oscillatory superfluid flow is determined by ^Jk/co (analogous to 
the viscous penetration depth of oscillating classical flow). This is so for various reasons. 
Firstly, the oscillation amplitude at the critical velocity v c ~ \/kco is a c = v c /a> ~ ^Jic/(D. 
Secondly, the intervortex spacing at v c is l c = l/\/L^ ~ sjKJG} (and hence scales as a c ). 
Moreover, and very interestingly, also from the dispersion relation of Kelvin waves, namely 
a>(k) = Kk 2 \n(2/kao) j\%, we obtain a wavelength X = 2n/k ~ \Jk/& when considering 
the logarithmic term to be approximately constant (w 12 in our case). It is also interesting 
to note that the group velocity v g ~ \/kco scales as v c . Therefore, it appears plausible that 
Kelvin waves traveling along a vortex that is attached to the oscillating sphere will be ex- 
cited very efficiently. On dimensional grounds we set a w °c a>~ 3 ' 2 °c (fc/co) 3 ' 2 to eliminate 
the time. It follows that another length scale, independent of CO, must be involved in order 
to get the dimensions right, probably the radius r of the sphere: a w ~ (K"/co) 3 ' 2 /r 2 . The 
average displaced liquid volume is AV « Ttr 2 ■ a w so that AV ~ (fc/a)) 3 ' 2 . Unfortunately, it 
remains unclear whether AV is, by any means, related to a 3 or Z 3 or A 3 ? 

There is a possible alternative way of reasoning based on diffusion of vorticity [ 1 2 1 . But 
it remains unclear as well, how in detail the intervortex spacing at the end of a turbulent 
phase increases from l c to the average value l c + a w by means of a diffusion process with a w 
to depend on the period T of the oscillation: a w °c 0)~ 3 ' 2 ex t 3 ' 2 . 



5 Conclusion 

We have a achieved detailed picture of the lifetime of the turbulent phases in terms of a fluc- 
tuating vortex density L. We have found that the fluctuations of L have a normal distribution 
and have determined the variance. The lifetimes of the turbulent states increase very rapidly 
with increasing drive but never diverge. Thus, turbulence is inherently unstable in our exper- 
iment. Our picture of the breakdown of potential flow by some remanent vorticity is not yet 
complete. While the occurrence of the Rayleigh distribution of the excess amplitudes Aa is 
plausible, the dependence of the rms value on the oscillation frequency cannot be explained 
yet. Also we do not have a sufficient set of data for a definite analysis of this effect. It would 
be very helpful to have a theory of oscillating superflow which up to date is not yet available. 
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